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Monte Carlo simulations are used to investigate the tricritical point properties of a 2d spin fluid. 
Measurements of the scaling operator distributions are employed in conjunction with a finite-size 
scaling analysis to locate the tricritical point and determine the directions of the relevant scaling 
fields and their associated tricritical exponents. The scaling operator distributions and exponents 
are shown to match quantitatively those of the 2d Blume-Capel model, confirming that both models 
belong to the same universality class. Mean-field calculations of the tricritical point properties are 
also compared with the simulation measurements. 
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I. INTRODUCTION 



For tricritical phenomena, the highest dimension in which non-classical behaviour can be observed is d — 2 [Q. 
Consequently, 2d tricritical phenomena has been the subject of a large number of previous investigations, employing 
a wide variety of techniques, including series expansions mean-field theory renormalisation group (RG) jfj-f), 
transfer matrix [p|-p"2"||, Monte Carlo simulations (P^,[l4|] and Monte Carlo RG methods jl5] [Dj. To date, however, 
this interest has focused almost exclusively on lattice-based spin models such as the Blume-Capel model or the spin- 
i next-nearest-neighbour Ising model. In this paper, we report the first detailed simulation study of 2d tricritical 
behaviour in an off-lattice spin fluid model. 

The model we consider is a simplified representation for a liquid of two-state molecules, and has been the subject 
of a number of previous studies in both its classical and quantum regimes ]isj] . In the present work, however, we shall 
consider only the classical limit, for which the configurational energy is given by: 

N N 

$({r, s}) = -J(r^) E s ^ + E U ^i) 

with Si = ±1 and where U(rij) is chosen to be a hard disk potential with diameter a. The distance-dependent spin 
coupling parameter J{rij) is assigned a square well form: 

J(r) = J a <r < 1.5a 

J(r) = elsewhere (1-2) 

The phase diagram of this model is expected to correspond to the situation depicted schematically in figure [l]. For 
high temperatures, there exists a line of Ising critical points (the so-called 'critical line') separating a ferromagnetic 
fluid phase from a paramagnetic fluid phase. The particle density varies continuously across this line. As one follows 
the critical line to lower temperatures, however, the size of the particle density fluctuations grows progressively. 
Precisely at the tricritical point, the fluctuations in both the particle density and magnetisation are simultaneously 
divergent. Lowering the temperature still further results in a phase separation between a low density paramagnetic 
gas and a high density ferromagnetic liquid. For subtricritical temperatures, the phase transition between these two 
phases is first order. 

Owing to the interplay between the density and magnetisation fluctuations, the tricritical properties of the spin 
fluid system are expected to differ qualitatively from those on the critical line. General universality arguments 
p9| predict that for a given spatial dimensionality, fluids with short-ranged interactions should exhibit the same 
tricritical properties as lattice-based spin systems. However, since fluids possess a continuous translational symmetry 
that lattice models do not, this proposal needs be checked. Additionally, experience with 'ordinary' (Ising) critical 
behaviour in simple fluids such as the Lennard- Jones fluid |2^j2^], shows that the reduced symmetry of fluids can 
profoundly influence certain non-universal aspects of the critical properties. Principal among these, are the directions 
of the relevant scaling fields associated with the fixed point, and the distribution functions of observables such as the 
particle density and energy. It is thus of interest to assess the extent of these 'field-mixing' effects in the tricritical 
fluid and to compare it with the situation at the liquid- vapour critical point of simple fluids. 

An accurate determination of the universal forms of the tricritical scaling operator distribution is also of considerable 
value. Such distributions are unique to a universality class and hence knowledge of their forms would be of considerable 
practical utility to future simulation studies of 2d tricriticality, serving as they do to simplify the computational task of 
locating the tricritical parameters. Moreover, as we shall see, the forms of the scaling operator distribution functions 
can impart important physical insight into the nature of the tricritical fluctuations. 

Our paper is broadly organised as follows. In section || we describe the finite-size scaling methods and other 



computational techniques employed in the study. We then proceed in section [II to detail the application of these 
techniques to Monte Carlo simulations of both the 2d spin fluid model described above, and the 2d Blume-Capel. The 
simulations yield accurate estimates of the location of the tricritical point for both models, as well as the universal 
forms of the tricritical scaling operator distributions and the directions of the relevant scaling fields. In the case of 
the spin fluid model, the estimates for the tricritical point parameters are compared with the results of a mean field 



calculation. Finally section IV details our conclusions. 



II. BACKGROUND 

The techniques we employ in this work have been previously developed in the context of simulation studies of Ising 
critical phenomena in a variety of fluid models, including a decorated lattice gas model [G3,p3|, a lattice model for 
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polymer mixtures p4j, and both the two and three-dimensional Lennard- Jones fluids pO| , pi| . In common with the 
strategy pursued in these previous works, we have chosen to work within the grand canonical ensemble, use of which 
affords effective treatment of the particle density fluctuations which are a central feature of fluid critical behaviour. 

Let us assume our system to be contained in a volume L d , with d = 2 in the simulations to be described later. The 
grand partition function is given by 



£ £11 { d A e-^ s})+ ^ N+H ^ ^ (2.1) 



N 

z L = ^^ 

N=0 { Si } i 



where N is the particle number, (3 = (KbT) -1 is the inverse temperature, p, is the chemical potential and H is the 
uniform applied magnetic field. 

The observables of chief concern to the present study are the (reduced) particle density 

p = L- d Na d , (2.2) 

the configurational energy density (which we express in units of J) 

u = L- d J- 1 mf,s}), (2.3) 

and the magnetisation. 

m. = L- d J2 s * ( 2 - 4 ) 

i 

The coarse-grained behaviour of the system in the vicinity of the tricritical point is controlled by three relevant 
scaling fields |j],^5|,^6| which we denote g, A and hi . In general these scaling fields are each expected to comprise linear 
combinations of the three thermodynamic fields T, p, and H [^7],|l2| ■ For the spin fluid model considered in this work, 
however, the configurational energy is invariant with respect to sign reversal of the spin degrees of freedom. This 
special symmetry implies that the tricritical point lies in the symmetry plane H = 0, and that the scaling field h! 
coincides with the magnetic field H, being orthogonal to the p, — T plane containing the other two scaling fields, g 
and A. Thus we can write 

ti = H-H t (2.5a) 
A = (p - ih) + r(T - T t ) (2.5b) 
g = T-T t + s(p- fi t ) (2.5c) 

where the subscript t signifies tricritical values and the parameters s and r are system-specific 'field mixing' parameters 
that control the directions of the scaling fields in the fi-T plane. The scaling fields g and A are depicted schematically 
in figure |l|(b) . One sees that g is tangent to the coexistence curve at the tricritical point Q , so that the field mixing 
parameter r may be identified simply as the limiting tricritical gradient of the coexistence curve. The scaling field 
A, on the other hand, is permitted to take a general direction in the fi-T plane which does not necessarily have to 
coincide with any special direction of the phase diagram |27j] . 

Conjugate to each of the scaling fields are scaling operators, defined by the requirements 

(M) = L- d d\nZ L /dti (2.6a) 
(V) = L- d d\nZ L /d\ (2.6b) 
(£)=L- d d\nZ L /dg (2.6c) 



from which it follows (utilising equations 2.1 and 2.5a — 2.5c), that 



M = m (2.7a) 

£>=-—!— [p-su] (2.7b) 
1 — sr 

£ = —!—iu~rp} (2.7c) 
1 — sr 

Motivated by our experience with ordinary critical phenomena in simple fluids [ p0[ we make the following finite-size 
scaling ansatz p8[ for the limiting (large L) near-tricritical distribution of pl (p, u, m) 
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p L (p, u, m) ~ j^—pLia^L^M, a£ 1 L d - ya V, a^L^S, ai L Vl ti, a 2 L y2 \, a 3 L y3 g) 



(2. 



where Pl is a universal scaling function, the are non- universal metric factors and the yi are the standard tricritical 
eigenvalue exponents (^9|. P reci sely at the tricritical point, the tricritical scaling fields vanish identically and the last 
three arguments of equation 2.8 can be simply dropped, yielding 



PL(p,U, m) 



1 



1 — sr 



plia^L^M, a^L d - y2 V, a^L d ^E) 



(2.9) 



where p* L is a universal and scale invariant function characterising the tricritical fixed point. 

In what follows we shall explicitly test the proposed universality of equation 2.9 for the case of the spin fluid model, 
by obtaining the form of p* L and comparing it with that for the tricritical 2d Blume-Capel model. 



III. RESULTS 



A. computational aspects 

The Monte-Carlo simulations of the spin fluid model were performed using a Metropolis algorithm within the grand 
canonical ensemble. Particle insertions and deletions were carried out using the prescription of Adams |^[n]. When 
attempting particle insertions the spin of the candidate particle was randomly assigned the value +1 or —1 with equal 
probability. Spin flip attempts were performed at the same frequency as insertion/deletion attempts. 

In order to facilitate efficient computation of interparticle interactions, the periodic simulation space of volume 
L 2 was partitioned int o I 2 cubic cells each of side 1.5a, corresponding to the interaction range of the interparticle 



potential (cf. equation IT). This strategy ensures that interactions emanating from particles in a given cell extend 
at most to particles in the 8 neighbouring cells. We chose to study three system sizes corresponding to I = 12, 16 and 
20, containing, at coexistence, average particle numbers of (N) = 120, 210 and 330 respectively. For the I = 12, 16 
and 20 system sizes, equilibration periods of 10 5 Monte Carlo transfer attempts per cell (MCS) were utilised, while 
for the I = 16 and / = 20 system sizes up to 2 x 10 6 MCS were employed. Sampling frequencies ranged from 25 MCS 
for the I — 12 system to 100 MCS for the I = 20 system. Production runs amounted to 1 x 10 7 MCS for the I = 12 
and up to 5 x 10 7 MCS for the I = 20 system size. At coexistence the average acceptance rate for particle transfers 
was approximately 16%, while for spin flip attempts the acceptance rate was approximately 5%. 

During the production runs, the joint probability distribution pl(p, U, m) was obtained in the form of a histogram. 
In order to increase computational efficiency, the histogram extrapolation technique was employed. Use of this 
technique permits histograms obtained at one set of model parameters to be reweighted to yield estimates appropriate 
to another set of model parameters. The method is particularly effective close to a critical point where, owing to the 
large fluctuations, a single simulation permits extrapolation over the entire critical region. 

As an aid to locating the tricritical point, the cumulant intersection method was employed f33f . The fourth order 
cumulant ratio Ul is a quantity that characterises the form of a distribution | p4| . It is defined in terms of the fourth 
and second moments of a given distribution 

Ul = 1- < > (3.1) 
3 < m z > z 



The tricritical scale invariance of the distributions pl(T>),pl(A4) and Pl(S), (as expressed by equation 2.9), implies 
that at the tricritical point (and modulo corrections to scaling), the cumulant values for all system sizes should be 
equal. The tricritical parameters can thus be found by measuring Ul for a number of temperatures and system sizes 
along the first order line, according to the prescription given below. Precisely at the tricritical temperature the curves 
of Ul corresponding to the various system sizes are expected to intersect one another at a single common point. 



B. The 2d Blume-Capel model 



In seeking to confirm the proposed universality linking the tricritical point of the 2d spin fluid model to those of 
2d lattice models, it is first necessary to determine the tricritical operator distribution functions for a simple lattice 
model. To this end we have also performed a simulation study of the Blume-Capel model on a periodic square lattice, 
the Hamiltonian of which is given by 
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h = -j J2 s iSj +Dj2si+ H J2 si ( 3 - 2 ) 



<l,J> 



with Si = —1,0, 1. Here H is a uniform magnetic field and D is the so-called 'crystal field'. As with the 2d spin fluid 
model, the symmetry of the configurational energy under spin sign reversal implies that the tricritical point lies in 
the symmetry plane H = 0. 



Previous MCRG 15|lq] investigations place the tricritical point of the 2d Blume-Capel model at KsTt/J = 
0.609(3), A = 1.965(15), while a more recent transfer matrix study §] gives K B T t /J = 0.610(5), A = 1.965(5). 
Using these estimates as an initial guide, we performed extensive Monte Carlo simulations of the model using a 
vectorised Metropolis algorithm on a Cray-YMP. Four system sizes of linear extent L — 12,20,32,40 were studied 
and following equilibration, runs ranging from 5 x 10 6 Monte Carlo sweeps (MCS) for the L—12 system, to 2 x 10 7 
MCS for the L=40 system were performed. The quantities measured in the course of these runs were 

p' = E s ?' ( 3 - 3a ) 

u = — J SiSj, (3.3b) 

<i,j> 

m = Si. (3.3c) 

Here we note that p' = X) s ? pl avs the same role as the density in the spin fluid model, being discontinuous across the 
first order line but continuous on the critical line. This fact is most clearly evident in the lattice gas representation 
of the Blume-Capel model, where the crystal field D appears as a chemical potential [jj5|. 

During the simulations, the joint distribution pl(p' ,u,m) was collected in the form of a histogram. To determine 
the locus of the first order line (and hence locate the tricritical point in which it formally terminates), the coexistence 
symmetry criterion for the operator distribution pl(D) was utilised. This criterion is the analogue for asymmetric 
first order transitions, of the order-parameter distribution symmetry condition applicable to symmetric first order 
transitions such as that of the subcritical simple Ising model [|23|. For a given temperature T, the first order transition 
point can thus be located by tuning the crystal field D and the value of the field mixing parameter s, within the 
histogram reweighting scheme, until the operator distribution pjjT>) is symmetric in T> — (T>). 

The first order line and its finite-size analytical extension EQ] was determined in this way for temperatures in 
the range KbT/J = 0.59-0.625, and for each of the 4 system sizes. The corresponding values of the cumulant 
ratio U"£ along this coexistence line are presented in figure | as a function of the temperature. To within numerical 
uncertainties the cumulant values for each system size intersect at a common temperature, which we estimate as 
KbTi/J = 0.608(1). The corresponding estimate for the tricritical field is D t = 1.9665(3). Clearly these values are 
in excellent agreement with the aforementioned estimates of previous studies. 

In the following subsection we compare the measured forms of the tricritical operator distributions Pl{-M),Pl('D) 
and Pl{£) of the 2d Blume-Capel, with those of the 2d spin fluid model. 



C. The spin fluid model 



The procedure for locating the tricritical point of the 2d spin fluid model followed the same pattern as that for the 
2d Blume-Capel model described above, except that in the present case no prior estimates for the tricritical point 
were available. It was thus necessary to search for the approximate location of the tricritical point by performing 
a number of short runs in which a temperature was chosen and the chemical potential tuned. Observations of 
the behaviour of the density from these short runs suggested that the tricritical point lay close to the parameters 
K B T/J = 0.58, j3fi = -1.915. 

Having obtained an approximate estimate of the location of the tricritical point, long runs were carried out for each 
of the three system sizes I — 12, 16, 20. As with the Blume-Capel model, the coexistence symmetry condition was then 
applied to the operator distribution pl (T>) in conjunction with the histogram reweighting scheme, in order to determine 
the first order line and its analytic extension in the fi—T plane. The locus of this line is shown in figure while the 
measured values of Jf£ along this coexistence line are shown in figure [| for the three system sizes I — 12, 16, 20. 
Clearly the curves of figure [| have a single well-defined intersection point, from which we estimate the tricritical 
temperature as being KsT t /J = 0.581(1). The associated estimate for the chemical potential is (3pt — —1.916(2). 
Typical near-tricritical configurations for the I = 20 system are shown in figure 

In figure ^| we present the forms of the operator distributions Pl{M), pl(2?) and Pl{£) corresponding to the 
designated values of the tricritical parameters. The value of the field mixing parameter r implicit in the definition 
of £, was assigned the value r — —2.82, as obtained from the measured gradient of the phase boundary in the fi-T 
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plane at the tricritical point. The value of the field mixing parameter s was assigned, as previously described, so that 
Pl(T>) satisfied the symmetry condition. However, the resulting estimates of s were found to exhibit a systematic 
finite-size dependence. This effect has also been previously noted (albeit with much reduced magnitude) in a recent 
study of critical phenomena in the Lennard- Jones fluid, and is traceable to the finite-size dependence of the average 
critical energy p3|po[ . For the three system sizes, I = 12, 16, 20, we found s = —0.031, —0.020, —0.013 respectively. 
Interestingly, these values are at least an order of magnitude smaller than those measured at the critical point of 
the 2d and 3d Lennard- Jones fluid, a finding which we discuss further in section [V. This smallncss implies that the 
scaling field A almost coincides with the p axis of the phase diagram. 

Also included in figure ^| are the measured tricritical operator distributions for the 2d Blume-Capel model. In 
accordance with convention, all the operator distributions have been scaled to unit norm and variance. Clearly in 
each instance and for each system size, the operator distributions collapse extremely well onto one another as well as 
onto those of the tricritical Blume-Capel model. 

The measured scaling operator distributions also serve to furnish estimates of the eigenvalue exponents j/i , 1/2 and 
U3 characterising the three relevant scaling fields. These exponents are accessible via the respective finite-size sca ling 
behaviour of tricritical distributions of pl{M), PlCD) and pl{£)- Specifically, consideration of the scaling form 2.9 
shows that the typical size of the fluctuations in a given operator O vary with system size like SO ~ L~( d ~ Vi K 
Comparison of the standard deviation of a given operator distributions as a function of system size thus affords 
estimates of the appropriate exponents. From the measured variance of the spin fluid operator distributions, we find 
Ui = 1.93(1), 2/2 = 1-80(1), 2/3 = 1.03(7). As far as 2/1 and y 2 are concerned, these estimates are in excellent agreement 
with exact conjectures 2/1 = 77/40, 2/2 = 9/5 [p6|-|3^]. The situation regarding the eigenvalue exponent 2/3, on the other 
hand, is not so satisfactory, there being a sizable discrepancy with the exact value of 2/3 = 4/5. This discrepancy 
stems, we believe, from two sources. Firstly, since the operator distribution Pl{£) is highly sensitive with respect 
to the designation of the value of the field mixing parameter r implicit in the definition of £, small uncertainties in 
the estimate of r can lead to significant larger errors in the measured variance of pl{£)- A similar effect was also 
previously observed at the liquid- vapour critical point of the Lennard- Jones fluid |||] . Secondly, and as we discuss in 
section IV, the near Gaussian character of Pl(£) signifies the absence of strong fluctuations in the £, in which case it 
is questionable whether a finite-size scaling can be reliable applied to Pl(£) at all. 

In view of this problem we have adopted a rather different approach for measuring 2/3 based on the scaling properties 
of U]?, a quantity which does show strong fluctuations and which is also insensitive to the designation of the field 
mixing parameter r. It can be shown |33|, that the maximum slope of the cumulant ratio near T c varies with 
system size like L Va . Using the histogram extrapolation technique, we have obtained the temperature dependence of 
this slope for the spin fluid model. The results yield the estimate 2/3 = 0.83(5), which agrees to within error with the 
exact conjecture. 

Turning now to the observables m, p and u, we plot in figure the measured distributions of these quantities 
at the assigned tricritical parameters. Here we note that as is the case with ordinary critical phenomena in the 
Lennard- Jones fluid plf , the energy distribution pl(u) differs qualitatively in form from the operator distribution 
Pl(£)- This finding, the origin of which is explained in detail in reference |23"| , reflects the coupling of the tricritical 
energy fluctuations to the density, the latter of which are stronger and thus dominate for large L. The influence of 
this coupling is also discernible as a small asymmetry in the tricritical density distribution. For the average tricritical 
density we find pt = 0.374(1), while for the average tricritical energy density we find u t = 0.778(2). The average 
magnetisation is of course strictly zero on symmetry grounds. 

Finally in this subsection, table || summarises the measured values of the fourth order cumulant ratios U£ for each 
of the three scaling operator distributions at the assigned values of the tricritical parameters. 



1. Mean-field calculations 



In this section we describe a simple mean-field calculation of the tricritical parameters of the 2d spin fluid. 

Let p{p, p; h) denote the pressure in the system for a given chemical potential p and volume V in the thermodynamic 
limit, where the equilibrium density p is given by the ratio of the average number of particles < N > and the volume 
V, p = lim y ^ 00 (< N > /V). Then 

p(p, p; h) = min p/ [-f(p'; h) + pp'\ (3.4) 
where /(//; h) is the free energy per volume, 

f(p';h) = lim -±lntr [exp (-0H N (h))] . (3.5) 
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In the mean field approximation we assume an interaction between a spin si and an effective field (g/N) X)j>i s i = <Z TO > 
where q is the effective coordination number and the ./V-particle Hamiltonian is written as 

N 

HZf = E u ( r *j) - E IMQ + h] Sl (3.6) 

(i<j) i=l 

and the mean field free energy fhiF is 

/mf(p'; ft) - Jim Infr [exp (-()H& F {h))] (3.7) 



/mf(/o'; ft) = / c z(p') + Tnin-n 



qp'm 2 



In 2 cosh (/3(qm + h)) 



(3.8) 



fd{p') is the free energy of a classical system with Hamiltonian iJ c ; = X)(i<j) ^( r y ) an( i the second term on the right 
hand side of Eq. (3.8) reaches its minimum at 

m(ft) = tanh \j3 (qm(h) + h)] (3.9) 

Since the coordination number q in the fluid is not fixed we approximate the effective field on one particle by 

qm = nip 1 I d 2 rJ(r)g(r) (3.10) 



where the fluid correlation function g(r) is taken from the Percus Yevick solution for hard discs, which can be found 
numerically. 

Two phase coexistence between a gas phase at density p g and a liquid phase phase at density pi is given by equal 
pressure in the two phases 



P(p g ,^h) =p(p h ij,;h) 



(3.11) 



This condition determines the mean-field chemical potential \imf for phase coexistence via Eq. (3.4). Implementing 



this criterion numerically, we find T t 1 = 1.006 (3p t — —1.319. Clearly this estimate for T t seriously overestimates 
the measured tricritical temperature (T t = 0.581) showing that mean field calculations of this type cannot be relied 
upon to provide accurate tricritical data, at least for 2d systems. 



IV. CONCLUSIONS 



In summary we have demonstrated that the tricritical ordering operator distributions of the 2d spin fluid can be 
mapped into excellent correspondence with those of the 2d Blume-Capel model. The existence of such a mapping 
represents perhaps the most stringent test of universality. There can thus be little doubt that despite their very 
different microscopic character, the two systems do indeed share a common fixed point. 

With regard to the scaling operator distribution themselves, we note that the form of p* L {£) is (to within the 
precision of our measurements) essentially Gaussian, as evidenced by the very small value of the cumulant ratio 
= 0.003(3). This Guassian behaviour implies that the tricritical fluctuations in £ are extremely weak and is thus 
a consequence of the central limit theorem. This weakness is further manifest in the very small measured values of 
the field mixing parameter s, as well as in the near absence of asymmetry in the tricritical density distributions, a 
situation which contrast markedly with that of the 2d and 3d Lennard-Jones fluids, where much stronger field mixing 
effects are observed in the density distribution |pl]| . 

Finally, the tricritical form of p* L {M.) merits special comment. We note that the three-peaked form of this distribution 
differs radically from the universal magnetisation distribution of the critical Ising model, which is strongly double- 
leaked in two dimensions p3|p9[ . The existence of a three peaked structure for tricritical phenomena reflects the 
additional coupling that arises between the magnetisation and the density fluctuations. Specifically, the central peak 
corresponds to fluctuation to small density, which are concomitant with an overall reduction in the magnitude of the 
magnetisation (cf. figure ||). Were one, however, to depart from the tricritical point along the critical line, these 
density fluctuations would gradually die out and a crossover to a magnetisation distribution having the double-peaked 
Ising form would be expected. In future work we intend to investigate the nature of this crossover in detail. 
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TABLE I. The fourth order cumulant ratio for the tricritical fixed point operator distributions. 




0.348(3) 
0.574(2) 
0.003(3) 



FIG. 1. (a) Schematic phase diagram of the spin fluid in the T-p plane, (b) Schematic phase diagram in the p-T plane 
showing the directions of the relevant scaling field g and A. 



FIG. 2. The measured cumulant ratio U L for the 2d Blume-Capel model along the first order line and its analytic extension, 
determined according to the procedure described in the text. 

FIG. 3. The line of first order transitions and its finite-size analytical extension in the fi-T plane, obtained according to the 
procedure described in the text. 



FIG. 4. The measured cumulant ratio U L for the 2d spin fluid model along the first order line and its analytic extension 
determined according to the procedure described in the text. 

FIG. 5. Typical particle/spin configurations of the I — 20 spin fluid near tricriticality. Spins values of +1 are denoted by 
filled circles, and spin values —1 by unfilled circles. 



FIG. 6. The scaling operator distributions for the 2d spin fluid at the designated tricritical parameters for each of the three 
system sizes I = 12,16,20. (a) p\{M), (b), p* L (T)) (c) p* L (E). Also shown for comparison are the corresponding distribution 
measured for the tricritical L — 40 2d Blume-Capel model. All distributions are expressed in terms of the scaling variable 
a^ 1 L d ~ Vi (O — O c ) and are scaled to unit norm and variance. Statistical errors do not exceed the symbol sizes. 
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FIG. 7. The measured tricritical point forms of: (a). The magnetisation distribution pz,(m) (the data points have been 
suppressed for clarity), (b) The density distribution Pl(p)- (c) The energy density distribution pl(u). 
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